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We use gravitational lensing of fast radio bursts (FRB) to try to constrain primordial black 
hole dark matter with the latest CHIME catalog and make a forecast for future results. 

¢ Using the intrinsic burst width and a calibrated flux-ratio threshold per FRB. 

¢ Taking into account the uncertainty in dispersion measure and redshift relation. 

¢ We outline an algorithm to detect lensed FRBs and to validate the possibility of detection. 


On ® 


¢ We demonstrate how stacking repeating FRBs can improve the 
constraints, especially for lower masses. flux 

¢ We Identified one FRB with a double peak as a candidate for 
strong gravitational lensing. 


Introduction 


FRBs are bursts of millisecond duration in radio frequencies. 
As the light interacts with charged particles as it propagates 
through the universe, the more energetic photons can push 
through the free electrons with little effect on their speed. 
This results in dispersion of the burst as seen in the figure to 
the Right [1]. Since the dispersion measure (DM) is 
proportional to the distance traveled, it can be used to 
estimate that distance and the redshift of the source z;. We PLIES SESS SOB eth ser ios 
use the Macquart relation [2] to estimate Zz; from the dispersion | en i er! 
measure while taking into account the uncertainty in their relation. pty aranaaee tal) 

We model the DM as consisting of: 


Where DMyy is the total Milky-Way contribution and the extra-galactic DM contribution is 


DMpost + DM 
DMeg = DM cosmic + ts 
DM-osmic is the contribution of the intergalactic medium (IGM), and DM,,., and DM,,- 


result from the FRB host galaxy and source environment. 


Observed frequency (MHz) 


The probability that an FRB is lensed is approximately the lensing optical depth. 


Zs 


T(M_, Zs) = | dy(z,)(1 + z,)*n,o(M,, Z1) 


0 
Where y(z) is the comoving distance at redshift z and n,; is the comoving number density of 
the lens. a(M,, Z,) is the effective cross section for lensing, it is an anulus defined by two 
requirements: 
¢ The first is that the two images are completely separated, for this we require that the 
separation is greater than the width of the FRB. This gives us the minimal normalized 
impact parameter for lensing y,,;. 
¢ The second is that the second image is strong enough to detect, this gives a maximal flux 
ratio between the two images R,, we find it through our algorithm by simulating the 
bursts as lensed and finding the maximal R- that enables detecting the lensing, Re. From 
this we get the maximal normalized impact parameter for lensing y+. 


Summing over every source’s lensing optical depth gives the total optical depth of a lens with 
mass M,, T(M,). 
7M.) =) t(25 i. Mi, Wir Ry) 
L 
The sum is over all FRBs that pass our validation process and their optical depths are a 
function of their redshift, lens mass, intrinsic width and Re that we find for each FRB. 


If none of the events are lensed the bound we can place on the fraction fp, of dark matter 
allowed in the form of PBHs is fpy < 1/T. 


Repeaters 


For repeating FRBs we can stack their N repetitions, so 
that the signal gets magnified by a factor of N, while the os 
noise only gets magnified by a factor of VN, therefore 
allowing to increase Rreby a factor of VN and improve = 
the sensitivity to strong lensing by PBHs. This relieson — 
two assumptions: (i) one and all FRB repetitions are 
lensed by the intervening PBH; (ii) a meaningful fraction 
of the repetitions have small intrinsic widths so that they 
can be efficiently stacked together. 

This figure shows fp, for: 1000 non-repeating events (blue), 100 events that repeat 10 times 
(orange). We simplified and assumed all FRBs can be detected as lensed with Ry = 9, in 
reality this is not the case. The full details of these calculations can be found in our paper. 
The dashed line is 30M_, as you can see, fewer repeating events can constrain fp, in low 
masses where many non-repeating events fail to give a constraint. 
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[1] Image from: Essential Radio Astronomy James J. Condon and Scott M. Ransom 
[2] Monthly Notices of the Royal Astronomical Society, Volume 509, Issue 4, February 2022, Pages 4775-4802 
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Source FRB 
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Before ruling out FRBs as lensed, we want to see for what Re 
they can be detected es lensed if they can indeed be detected. 
To do this, we simulate the events as lensed with different R¢ 
and find the maximal one that allows detection. 
The lensing detection is done through autocorrelation. 
First we define a threshold based on the autocorrelation of 
unlensed burst. For the unlesed burst we take the burst 
(highlighted in orange in the Fig. 1) and its preceding 
section of noise where we can be sure there is no echo from : > 40 = 
lensing. We then flip it to place the burst in the beginning Atims] 
(Fig. 2) and calculate its autocorrelation. Image 3 shows this autocorrelation. The beginning 
of the autocorrelation is high because of the burst overlapping with itself, we discard this 
part and only use the rest (the shaded part of Fig. 3) to calculate the mean pu and the 
standard deviation o of the autocorrelation. We set the threshold at uw + 2.3250 
(corresponding to a detection of an outliner at 99% — C.L). 
To simulate the lensing add the echo as follows: we subtract the mean of the noise from the 
burst, divide it by Rr and shift it forward in time enough to ensure a full separation between 
the burst and the echo. The lensed signal is shown in Fig. 4 with the echo highlighted in 
orange. Fig. 5 shows the autocorrelations of the FRB in three cases: lensed with Rr = 5, 
lensed with Ry = 10, and unlensed. For an FRB to be considered valid at a certain Rr we 
require the autocorrelation to pass the threshold (blue line in Fig.5) in the place we added 
the echo. 
An FRB is considered as a candidate for lensing if the autocorrelation of the original signal 
passes the threshold past the dashed line. This line marks the beginning of the lensing 
search, to its left the autocorrelation is high due to the burst overlapping with itself. 
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Results 
We focus on the range 10 — 100M o@ as it is motivated by oe ; “ 
LIGO and as there already are sufficient constraints above ,,... - * ee # °° 7 
100M_. | SAS. aay, . : 
We first make two gross cuts to the CHIME FRB data. offo” ates 


Width[s] 


We filter out FRBs with widths => 6ms which are too 
wide to constrain PBHs in the desired mass range, and 
FRBs with SNR < 25 as the noise in these FRBs is simply 
too large to allow a detection of a lensed echo. 10 
This cut leaves us with only 143 FRBs to be validated using the SNR 

process described above. The 114 FRBs that pass our validation process are shown in the 
figure above with their width, SNR and Re. 


Lensing candidate 


We found two candidates for lensing  *| 
and were able to rule one out due to ; 
different peak structure. The second 
candidate is FRB20190627B, here are 
Its light curve and autocorrelation. 
This double burst could be explained 
via strong lensing bya O(10M._) PBH. » 
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Forecast 


With 1 FRB at most being lensed the 114 FRBs 

that pass the validation process are not enough 

to place a bound on fpy, as you can see in this 
figure, it results in fpy > 1 (red). , 
However with only three times the current data 3 
we can constrain fpy even at low masses 0.5 
(blue). The dashed line at 30Mo is shown for 
orientation. 
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